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Abstract 

A mathematical description of the refractoriness period in neuronal 
diffusion modeling is given and its moments are explicitly obtained in a 
form that is suitable for quantitative evaluations. Then, for the Wiener, 
Ornstein-Uhlenbeck and Feller neuronal models, an analysis of the fea- 
tures exhibited by the mean and variance of the first passage time and of 
refractoriness period is performed. 

1 Introduction 



Mathematical descriptions of dead time, or refractoriness, in neuronal modeling 
have long traditions dating back at least to mid sixties when special attention 
was devoted to the description of the evolution of networks of switching ele- 
ments whose behavior was meant to simulate that of physiological neurons via 
certain rather drastic simplifications 0|. Furthermore, the approach to neural 
modeling was shown to bear certain strong analogies with the stochastic de- 
scription of the input-output features of radioactive particle counters. In such 
context, as early as 1948, W. Feller proved that under a suitable formulation 
all problems concerning single counters reduce to special instances of the theory 
of the summation of random variables. Exploiting the above mentioned analo- 
gies, the simplest neural model may be conceived as a black box possessing the 
following distinctive features: (i) it is a threshold element, (ii) its output re- 
sponse consists of pulses of constant amplitude and width and (in) there exists 
a constant dead time. More accurately, one could define this dead time also as a 
deterministic function of certain measurable parameters, such as time or input 
pulse amplitude, or view it as a stochastic process. 

As a first attempt towards a quantitative treatment of the dead time ef- 
fects in neural modeling, we look at the input of the neuron as a randomly 
distributed Poisson-type pulse train. Its output is then determined by imposing 
the restriction that following each input pulse a dead time period is activated 
during which no further pulses can be produced at the output. Even for such 
an oversimplified instance the investigation of the role played by the dead time 
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in determining the distribution of the output when the input is described by a 
given distribution is a very challenging task. 

Let t denote this dead time, i.e. the time interval following every firing 
during which the neuron cannot fire again. Let us assume that the net input to 
the neuron in time interval (0, T) is modeled by a sequence of positive pulses of 
standard strength whose time of occurrences are Poisson distributed with rate 
A > 0. We purpose to determine the distribution I1„(T, r) of the output pulses 
as a function of dead time r. A rather cumbersome amount of computations 
leads one to conclude that the assumed input distribution 

(1) Pn(T) = e- XT , T>0, n = 0,1,2,... 

generates the following firing distribution valid for all n > 1 (cf. 0]) : 

U n (T,r) = &[T - (n - l)r] {l - g ~ ~ " 1" " ~ } 

n 

\ (T-nr) 



fc=0 

fc 

\ ( / . - X — > - * I 1 -nil 



X k (T - nr) 



fc! 

fc=0 



(2) -tf(T-nr) 
where denotes the Heaviside unit step function: 

( 3 ) = { J; * < 2. 

Although the stated problem has been the object of several investigations (see, 
for instance, [Hj), a quantitative evaluation of the effect of dead time on the 
statistical parameters of the output appears to be still lacking. 

In the remaining part of this paper, we shall outline a totally different ap- 
proach towards the inclusion of refractoriness in the neuronal model. As in 
P| and P], we model the time course of the membrane potential by a time- 
homogeneous diffusion process and then assume that the firing threshold acts 
as some kind of elastic boundary characterized by preassigned reflection and 
absorption parameters. In other words, we assume that an action potential is 
released whenever the process first attains the firing threshold. After the firing, 
a period of refractoriness of random duration occurs, at the end of which the 
process is instantaneously reset at a fixed state. Then, the subsequent evolution 
of the action potential proceeds as before, until the threshold is again reached. 
A new firing then occurs, followed by a new period of refractoriness, and so on. 
Use of the above approach allows one to mimic the effects of refractoriness for 
the specified neuronal model. 

In order to be able to apply the specified paradigm to the description of 
neuronal models in the presence of refractoriness, an investigation of certain 
general features of diffusion processes in the presence of an elastic boundary is 
necessary. This task will be accomplished in Section 2, where we shall analyze 
the features of the moments of the random variable modeling the neuron's in- 
trinsic refractoriness. In Section 3, a specific analysis will be provided of three 
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neuronal models based on the Wiener, Ornstcin-Uhlenbeck and Feller diffusion 
processes, and a comparative discussion of the refractoriness features exhibited 
by these models will be performed. 



2 Effect of Refractoriness 



Let {X(t),t > 0} be a regular, time-homogeneous diffusion process, defined over 
the interval / = (ri, r 2 ), characterized by drift and infinitesimal variance A\(x) 
and A 2 (x), respectively. Throughout, we shall assume that Feller conditions on 
these functions are fulfilled [3]. Let h(x) and k(x) denote scale function and 
speed density of X(t): 

h(x) = expi-2 f 1 dz\, k(x) — 

17 A 2 (z) y A 2 {x)h{x) 

and 

rv rv 
H(r u y}= h(z)dz, K(n,y]= k{z) dz 

Jr\ J T\ 

scale and speed measures, respectively. 

We define the random variable "first passage time" (FPT) of X(i) through 
S (S € I) with X(0) = x<S: 

(4) T x = inf{i : X(t) > S}, X(0) = x < S. 
Then, 

(5) g(S,t\x) = -^P(T<t), x<S 

is the FPT pdf of X(t) through S conditional upon X(0) = x. 

In the neuronal modeling context the state S represents the neuron's firing 
threshold and g(S,t\x) the firing pdf. 

More realistically then in past approaches, here we shall assume that after 
each firing a period of refractoriness of random duration occurs, during which 
either the neuron is completely unable to respond, or it only partially responds 
to the received stimulations. To this end, we look at the threshold S as an elastic 
barrier being 'partially transparent', in the sense that its behavior is intermedi- 
ate between total absorption and total reflection. The degree of elasticity of the 
boundary depends on the choice of two parameters, a (absorbing coefficient) 
and (3 (reflecting coefficient), with a > and j3 > 0. Hence, p a := (3/(a + (3) 
denotes the reflecting probability at the boundary S, and f — p R = a /(a + 0) 
the absorption probability at S. We denote by T x the random variable de- 
scribing the "first exit time" (FET) of X(t) through S if X(Q) = x < S, and 
by g e (S,t\x) its pdf. The random variable T r will denote the "refractoriness 
period" and g r (S,t\S) its pdf. Since T x can be viewed as the sum of random 
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Figure 1: Illustrating first passage time t x through S (i.e. firing time), neuron's 
refractoriness period 5 r and first exit time t x for a single sample path x(i) of X{t). By 
t x , t x and 8 r we have indicated the appropriate values of T x> T x and T r , respectively. 



variable T x describing the first passage time through S (firing time) and of T r 
(see Figure^ one has: 

(6) g e (S,t\x)= g{S,T\x)g r {S,t\S,T) dr. 

Jo 

In the sequel we assume that one of the following cases holds: 

(i) ri is a natural nonattracting boundary and K(ri, y] < +oo; 

(ii) ri is a reflecting boundary or it is an entrance boundary. 

Under such assumptions, if x < S the first passage probability P(S\x) from x 
to S is unity and the FPT moments t n (S\xo) = E(T") are finite and can be 
iteratively calculated as 

t n (S{2) := / t n g(S,t\x)dt = n h(z) dz k{u) t n -i(S\u) du (n = 1,2, . . .), 

JO Jx Jrx 

where ^o(S'|a;) = P(S\x) = 1 (cf., for instance, [7]). 

Theorem 2.1 Under the assumption (i) and (ii), if a > the first exit time 
probability 

r+oo 

P(S\x) := / g e (S,t\x) dt (x < S) 
Jo 

is unity. 

Proof: We consider separately the cases (i) and (ii). 

Case (i) Let P(Si,S\x) (ri < Si < x < S) be the first exit time probability 
through the elastic boundary S in the presence of an absorbing boundary Si. 
This is solution of the differential equation 

(8) d^ + M^d*Mx) =Q 

dx 2 dx 2 
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subject to conditions 

(9) lim ib (x) = 0, a lim[l - ip (x)] - & limjfe- 1 ^) , ) = 0- 
Since 

#o(a) , A 2 (x) d 2 i/) (x) 1 d \ 1 dip (x)~ 

(10) A x (af) 



dx 2 c&r 2 fc(ir) cfa: [M^O dx 

from JSJl one has 

(11) ipo(x)=A + Bj h(z)dz, 

where A and B are arbitrary real constants. By imposing boundary conditions 
©, one obtains 



(12) P(S u S\x) = 



a / h(z) dz 



s x 



s 

a I h{z) dz + (3 

Si 



Since r*i is a natural nonattracting boundary one has H(r\,x] — +oo; hence, 
making use of one has 

P(S\x) := lim P(5i,5|a;) = 1, 

where the last equality follows by exploiting l'Hospital's rule. 
Case (ii) Let P(S\x) (n < a; < 5*) the first exit time probability through 
the elastic boundary S in the presence of a reflecting boundary or of an en- 
trance boundary r\. This is solution of the differential equation JSJ subject to 
conditions: 
(13) 

l]m u-x {x) ^oMi = a lim[1 _ Mx)] _ p lim U-i {x) ^oM} = o. 

i|r,l dx ) xtS L xtSl dx J 

Since (|10() holds, from JSJ) one obtains again the general solution (|llf) . By 
imposing boundary conditions (|13|) . one finally has = 1. 

■ 

Theorem 2.2 Under the assumption (i) and (ii), if a > the first exit time 
moments t n (S\x) = E(T™) can be iteratively calculated as 

/■oc ( />5 pZ 

t n (S\x) := / i n e (Sailor) dt = n < / h(z)dz I k{u)t n -\{S\u) du 

JO \^Jx Jri 

(14) +- / k(u)t n -i(S\u) du\ (n = 1,2,.. .; x < S), 



where to(S\x) = P(S\x) = 1. 
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Proof: A derivation of (|14|l follows from the properties of elastic boundaries. 
We consider again separately the cases (i) and (ii). 

Case (i) Let t n (S\,S\x) (r± < Si < x < S) be the first exit time moments 
through the elastic boundary S in the presence of an absorbing boundary S\. 
This is solution of the differential equation 



(15) 



dip n (x) A 2 (x) d 2 tb n (x) 
Mx) —j— + — = = -n^-iW 



dx 



dx 2 



subject to conditions 

(16) lim^ n (a;) = 0, a YmiibJx) + (3 Yimih^ix) 1 = 0. 

xlSx x]S x^S I dx ) 

The general solution of l(T5|l is 

r-X pX pZ 

(17) tb n (x)=A + B / h(z) dz — n I h(z) dz / k(u) ip n -i(u) du, 



where A and B are arbitrary real constants. By imposing boundary conditions 
l|16[l. one has 



f„(5i,5|a;) 



a / h(z) dz + (3 



/ h(z) dz k(y)tn-i(Si,S\y) dy 



-a 



h(u) du ■ / h(z) dz / k{y) * n _i(Si, dy 



+ / h(u) du ■ I h(z) dz / k(y)t n _i(Si, S\y) dy 



(18) 



Since ri is a natural nonattracting boundary and -?T(ri,y] < +oo, making use 
of l(T^|l and by applying l'Hospital's rule, one has 



t n (S\x) := lim t n (Si,S\x) =n 





a 



fc(y) ? n _i(5|y) 



+ / h(z) dz k(y)t n -i(Si,S\y) dy + h(u) du ■ k(y)t n -i(Si, S\y) dy \ , 



that identifies with the the right-hand side of 1)14(1. 

Case (ii) Let t„(S'|x) (ri < x < S) be the first exit time moments through 
the elastic boundary S in the presence of a reflecting boundary or of an en- 
trance boundary r%. This is solution of the differential equation (|15f) subject to 
conditions 



(19) lim fe" 1 ^) = 0, 

xln dx 



a lim ib r,(x) + lim< h 1 (x) 
xl;S x\sl 



d4>n{x) ' 

dx 



= 0. 
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From l|15fl one obtains again the general solution l|17f) . By imposing boundary 
conditions (|19|l . one finally is led to (|14|) . 



Note that in the absence of refractoriness, (|14|) arc in agreement with J7J. In- 
deed, if j3 = one has t n (S\x) — t n (S\x). The following remark shows that 
FET moments t n (S\x) are related to the FPT moments t n (S\x). 

Remark 2.1 Under the assumption (i) and (ii), if a > one has 

(20) t n (S\x) = t n {S\x) W": 1 ) t n -x-j(S\x) [ & 'fc(«)?i6%) du, 

a j~^Q V J / Jri 

(21) t n (S\x)=t n {S\x)+n^^C n ~ r )tj(S\x) f k{u) t n - X -j{S\u) du. 

^ 7=0 ^ 1 

Proof: Making use of (|14|l . relations (|20|) and l|21|) immediately follow by in- 
duction. 



Setting n = 1 in (|20|l or in (|21|l . one can see that the mean of first exit time 
is given by 

(22) t 1 (S\x)=ti(S\x) + £ [ k(u)du (x < S). 

a J ri 

Furthermore, setting n = 2 in lj!2UI) and in (J2U), one can obtain two equivalent 
expressions for the second order moment of first exit time: 

P - ten \ f S 7 / N 7 o P rS 



t 2 (S\x) = t 2 (S\x) + 2 - h(S\x) / k(u)du + 2- / k(u) h(S\u) du 
a v J ri a J ri 

(23) =ta(5|a:)+2^t 1 (5|a5) / k{u)du + 2^ [ k(u)h{S\u) du. 



Hence, the variance V(S\x) of the first exit time is given by 
(24) V(S\x) =V(S\x)+ (- f k(u)du] +2^ f fc(u)ti(%) du, 



n 



where V(S\x) denotes the FPT variance. 

Theorem 2.3 Under the assumption (i) and (ii), if a > the refractoriness 
period is doomed to end with certainty and its moments can be iteratively cal- 
culated as 
(25) 

E(T?):= [ t n g r (S,t\S)dt = n^ [ fc(u) t„_i(S|u) du (n=l,2,...). 
Jo a Jri 
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Proof: Integrating both sides of (JBJ in (0, +00) one has 

" + OO 

g r (S,t\S)dt=l, 



implying that the refractoriness period is doomed to end with certainty. Fur- 
thermore, from © we also have: 

t n {S\x):= / t n g e (S,t\x) dt = / dt t" / g(S, t\x) g r (S, t\S, r) dr 
Jo Jo Jo 

+ OG p- (-00 

dr.g(5,T|x) / i n Sr (5, r) di 

/r 

r + 00 r + co 

= d T g(S,T\x) (T + d) n g r (S,$\S)d$ 

Jo Jo 

= E ("W,(%)£(^)- (n = l,2,...). 

Hence, 

(26) E(T^)=t n (S\x)-^2 ( )t n ^(S\x)E(T^) (n=l,2,...). 

3=0 V' 

We now proceed by induction. Setting n = 1 in l(2^|) one sees that E(T r ) = 
ti(S\x) — ti(S\x). Hence, on account of lt2"2l . (l2~5)l holds for n = 1. Furthermore, 
assuming that (|26|) hold for j — 1,2, ... ,n, the right-hand side of (|26|l for n + 1 

becomes: 

( n + 1 )t n+1 ^(S\x)E(T?) 

3=0 ^ 3 ' 

= t n+1 (S\x) - t n+1 (S\x) - (ra + 1) - 22 I .J i n _j-(S|a;) / fe(u) du 

a j=0 

(27^=(n + l)- / k(u) t„(S\u) du, 



where the last equality follows from (|20Jl . From (|26|l we note that the left-hand 
side of H27JI is equal to E(T 7 n+1 ). Hence, if (|25|) holds for an arbitrarily fixed n, 
it also holds for n + 1, which completes the proof. 



Comparing (|14|l and i|25|) we note that 



E{T?) EElim t n (S*|x). 
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In particular, from l|25l) the first two moments and the variance of the refrac- 
toriness period are seen to be: 



E(T r ) = - [ k(u) du 



(28) E{T?) = 2- f k(u) htSlu) du + 2 ( £ / jfe(u) du] 

a J ri \a J ri J 

V(T r ) = 2^ jf jfe(u) *i(S|u) du + ^ ^ fc(u)du) . 

Comparing the first and last of l|28|) with (|22(l and (|24|) . we have 

(29) ti(S\x)=h(S\x) + E(T r ), V(S\x) = V(S\x) + V(T r ), 

i.e. the mean (variance) of first exit time through S starting from x is the sum 
of the mean (variance) of first passage time through S starting from x and of 
the mean (variance) of the refractoriness period. 



3 Analysis of three neuronal models 

In order to embody some physiological features of real neurons, several alter- 
native models have been proposed in the literature (cf, for instance, jS], |H] 
and references therein). In this Section we shall investigate the behavior of the 
refractoriness period for the Wiener, Ornstein-Uhlenbeck (OU) and Feller neu- 
ronal models. We assume that all three neuronal models are restricted to the 
same diffusion interval I = [u, +oo), having set r*i = v. 





h(S\p) 


E(T r ) 




E(T r ) 


E(T r ) 


E(T r ) 


a 2 




V R = 0-1 




V R = 0-5 


V R = 0.9 


p R = 0.99 


10. 


3.073451 E+2 


6.294544 E4 


-2 


5.665090 E+3 


5.098581 E+4 


5.608439 E+5 


20. 


7.331871 E+l 


9.425701 E-f 


-0 


8.483131 E+l 


7.634818 E+2 


8.398300 E+3 


30. 


3.936016 E+l 


2.021650 E+ 


-0 


1.819485 E+l 


1.637537 E+2 


1.801290 E+3 


40. 


2.663797 E+l 


8.663807 E- 


-1 


7.797426 E+0 


7.017684 E+l 


7.719452 E+2 


50. 


2.007160 E+l 


4.966112 E- 


-1 


4.469501 E+0 


4.022551 E+l 


4.424806 E+2 


100. 


8.937578 E+0 


1.281821 E- 


-1 


1.153639 E+0 


1.038275 E+l 


1.142103 E+2 


200. 


4.225259 E+0 


4.617762 E- 


-2 


4.155986 E-l 


3.740387 E+0 


4.114426 E+l 


300. 


2.765483 E+0 


2.760995 E- 


-2 


2.484895 E-l 


2.236406 E+0 


2.460046 E+l 


400. 


2.055224 E+0 


1.961207 E- 


-2 


1.765086 E-l 


1.588577 E+0 


1.747435 E+l 


500. 


1.635207 E+0 


1.518666 E- 


-2 


1.366799 E-l 


1.230119 E+0 


1.353131 E+l 



Table 1: Wiener model with fi — —0.5 and a 2 — 10-i, 100-i (i = 1,2, . . . , 5), restricted 
to I — [u, +oo) with v = —80. In the second column we have listed the FPT mean 
t\{S\p) with S — —50 and g — —70. Instead, in columns three, four, five and six we 
have respectively listed the mean of refractoriness period for p R — 0.1, 0.5, 0.9, 0.99. 
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V(S\ P ) 




V(T r ) 




V(T r ) 




V(T r ) 


V{T r ) 




a 2 






Pn =0-1 




Pn =0.5 




Pn = 0-9 


p R = 0.99 




10. 


9.254218 EH 


-4 


7.681238 EH 


-5 


3.544044 E- 


-7 


2.629677 E+9 


3.148772 E- 


-11 


20. 


4.970295 Eh 


-3 


1.310444 EH 


-3 


1.819075 Eh 


-4 


6.818541 E+5 


7.161989 En 


-7 


30. 


1.390880 E- 


-3 


1.385660 EH 


-2 


1.541363 E- 


-3 


3.770806 E+4 


3.364468 E- 


-6 


10. 


6.265060 E- 


-2 


3.874901 EH 


-1 


4.027854 E- 


-2 


8.002658 E+3 


6.297560 E- 


-5 


50. 


3.519348 EH 


-2 


1.638408 EH 


-1 


1.652135 EH 


-2 


2.925226 E+3 


2.101676 EH 


-5 


100. 


6.821593 EH 


-1 


1.804893 EH 


-0 


1.742704 EH 


-1 


2.526670 E+2 


1.463751 EH 


-i 


200. 


1.506377 Eh 


-1 


3.008332 E- 


-1 


2.861029 E- 


-0 


3.818526 E+l 


1.958992 E+3 


300. 


6.426684 Eh 


-0 


1.168741 E- 


-1 


1.106754 EH 


-0 


1.440657 E+l 


7.086384 E+2 


400. 


3.542131 En 


-0 


6.147195 E- 


-2 


5.809411 E- 


-1 


7.471651 E+0 


3.597818 EH 


-2 


500. 


2.239493 En 


-0 


3.778974 E- 


-2 


3.567134 E- 


-1 


4.555481 E+0 


2.165615 En 


-2 



Table 2: For the Wiener model and for the same choices of parameters of Tabled 
in the second column we have listed the FPT variance V(5|g) with S = —50 and 
g = —70, whereas in columns three, four, five and six we have respectively listed the 
variance of refractoriness period for p a — 0.1, 0.5, 0.9, 0.99. 



3.1 Wiener model 

The Wiener neuronal model is defined as the diffusion process X(t) characterized 
by the following drift and infinitesimal variance: 



(30) 



A 2 (x) = a 2 



{ll e R,cr > 0), 



restricted to I — [is, +oo), where on the regular boundary x — v a reflecting 
condition is imposed. For such process scale and speed functions are 



%)=exp{-^}, fc(x) = -Jexp{^} 



2 jJL X ~ 



Furthermore, the mean of first passage time is 
(31) 

( {S-x){S + x-2v) 



h(S\x) = { 



S — x 



a 



f 2fi(S-v)\ r -Afx{x-v)\ 



2 ix (x — v) 



H = 



For the Wiener model (JSOJ with (x = -0.5, a 2 = 10 ■ i, 100 • i (i = 1, 2, . . . , 5), 
restricted to I = [u, +oo) with r\=v= —80, in the second column of Tabled 
and of Table[3we have respectively listed the mean t\{S\o) and variance V(S\g), 
numerically obtained via J7J) with S = -50 and g = —70. Note that the FPT 
mean and variance decrease with a 2 . Being j3/a — p a /(l — P R ), in Tabled 
and in Table |21 we have respectively listed the values of mean and variance of 
refractoriness period, numerically obtained via (|28() for p R — 0.1, 0.5, 0.9, 0.99. 
We observe that E(T r ) and V(T r ) increase with p R for any fixed a 2 . 
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h(S\ P ) 


E(T r ) 




E(T r ) 


E(T r ) 


E(T r ) 


a 2 




P R = o.i 




p R = 0.5 


p R = 0.9 


p n = 0.99 


10. 


9.862135 E+3 


9.901436 E- 


-41 


8.911293 E+42 


8.020163 E+43 


8.822180 E+44 


20. 


2.600359 E+2 


3.452097 E+20 


3.106887 E+21 


2.796199 E+22 


3.075818 E+23 


30. 


7.956655 E+l 


2.140293 E- 


-13 


1.926264 E+14 


1.733637 E+15 


1.907001 E+16 


10. 


4.273886 E+l 


4.978530 E- 


-9 


4.480677 E+10 


4.032609 E+ll 


4.435870 E+12 


50. 


2.853092 E+l 


3.149993 E-\ 


-7 


2.834994 E+8 


2.551494 E+9 


2.806644 E+10 


100. 


1.038152 E+l 


1.006196 EH 


-3 


9.055763 E+3 


8.150187 E+4 


8.965206 E+5 


200. 


4.525217 E+0 


4.073683 E- 


-0 


3.666314 E+l 


3.299683 E+2 


3.629651 E+3 


300. 


2.890905 E+0 


5.465816 E- 


-1 


4.919234 E+0 


4.427311 E+l 


4.870042 E+2 


400. 


2.123662 E+0 


1.839895 E- 


-1 


1.655905 E+0 


1.490315 E+l 


1.639346 E+2 


500. 


1.678216 E+0 


9.103197 E- 


-2 


8.192877 E-l 


7.373589 E+0 


8.110948 E+l 



Table 3: OU model with ■& = 5, q = -70 and a 2 = 10 ■ i, 100 ■ i (i = 1, 2, . . . , 5), 
restricted to I — [v, +oo) with v = —80. In the second column we have listed the 
FPT mean t\(S\p) with S = —50, whereas in columns three, four, five and six we have 
respectively listed the mean of refractoriness period for p R — 0.1, 0.5, 0.9, 0.99. 



3.2 OU model 

The OU neuronal model is defined as the diffusion process X(t) characterized 
by the following drift and infinitesimal variance: 

(32) A x [x) = -- (x- g), A 2 = a 2 (g e R, a > 0, <& > 0), 
v 

restricted to I — [u, +oo), where on the regular boundary x — v a reflecting 
condition is imposed. For such process the scale and speed functions are 



Mx)=exp{-^-^}, /c(z) = ^ cxp{- 



da 2 



2 gx ' 



Furthermore, the mean of first passage time is: 

2 fc \fS-Q 



ti{S\x) = d 



+ 00 

^(fc + l) (2fe + l)!! 



(33) 



-2-9 exp 



+ 00 

E 

k=0 



(iz-e) 3 ! ^ 2 fe 



2fc+2 



f \y- Q) ) 

X o 2 $ J 



^ (2fc + l)!! 



v 



x — g 

2k 

g 



2k+2 



1 



(2k + l)fc! 



S — g 



2k+l 



, 2fe+l 

Q 



For the OU model ^ with i? = 5, g = -70, a 2 = 10 • i, 100 ■ i (i = 1, 2, . . . , 5), 
restricted to I = [u, +oo) with v = —80, in the second column of Table and 
of Table 21 we have respectively listed the mean ti(S\g) and variance V(S\g), 
numerically obtained via Q with S = —50. Furthermore, in Tableland in 
Table 0] we have also listed the values of mean and variance of refractoriness 
period for p a = 0.1, 0.5, 0.9, 0.99. Similarly to the case of the Wiener model, 
for the OU model the FPT mean and variance decrease with a 2 ; furthermore, 
E(T r ) and V(T r ) increase with p R for any fixed a 2 . 
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V(S\ P ) 




V{T r ) 




V(T r ) 


V(T r ) 


V{T r ) 


a 2 






p B =0.1 




p R = 0.5 


P R = 0.9 


p R = 0.99 


10. 


9.713857 EH 


-7 


9.803844 E- 


-83 


7.941114 E+85 


6.432302 E+87 


7.783086 E+89 


20. 


6.554937 Eh 


-4 


1.191697 Eh 


-11 


9.652749 E+42 


7.818727 E+44 


9.460659 E+46 


30. 


5.898427 E- 


-3 


4.580854 E- 


-26 


3.710492 E+28 


3.005498 E+30 


3.636653 E+32 


10. 


1.654790 E- 


-3 


2.478576 E- 


-19 


2.007647 E+21 


1.626194 E+23 


1.967694 E+25 


50. 


7.239524 EH 


-2 


9.922473 E+14 


8.037191 E+16 


6.510124 E+18 


7.877250 E+20 


100. 


9.240940 EH 


-1 


1.029849 En 


-6 


8.216362 E+7 


6.643966 E+9 


8.037646 E+ll 


200. 


1.728177 E- 


-1 


4.589508 E- 


-1 


1.607888 E+3 


1.112524 E+5 


1.320047 E+7 


300. 


7.020488 Ei 


-0 


2.765114 E+0 


4.639613 E+l 


2.159884 E+3 


2.393706 E+5 


400. 


3.780330 EH 


-0 


6.379705 E- 


-1 


8.179088 E+0 


2.710374 E+2 


2.741284 E+4 


500. 


2.357829 En 


-0 


2.431146 E- 


-1 


2.784683 E+0 


7.339088 E+l 


6.787980 E+3 



Table 4: For the OU model and for the same choices of parameters of Table [IJ in the 
second column we have listed the FPT variance V(S\q) with S = —50 and g = —70, 
whereas in columns three, four, five and six we have respectively listed the variance of 
refractoriness period for p R = 0.1, 0.5, 0.9, 0.99. 

3.3 Feller model 

The Feller neuronal model is defined as the diffusion process X(t) characterized 

by the following drift and infinitesimal variance: 

(34) 

A\(x) = — - (x — q), A 2 {x) = 2£(x-u) (q, v £ R, g> v,d > 0,£ > 0). 
v 

defined in / = [v, +oo), where x = v is regular if g — v < and entrance if 
g — v > £ fl, whereas the boundary +oo is natural. For such process the scale 
and speed functions are: 

*(*) = exp{-^} ( X -,)-^^\ k( X ) = i exp{--^} ( X -^-^- 

The mean of the firing time can be calculated; for x < S one obtains 
(35) 

»™ - ^ h + |Q) t (s -^;r^ {n(s^')}' 

For the Feller model (JUJ with i9 = 5, g = -70, v = -80, ^ = 0.5 • % (i = 
1,2, ... , 10), in the second column of Tableland of TablcHJwe have respectively 
listed the mean t±(S\g) and variance V(S\g), numerically obtained via Q with 
S = —50. Note that the FPT mean and variance decrease with £. Furthermore, 
in Table [S] and in Table El we have listed the values of mean and variance of 
refractoriness period for p R = 0.1, 0.5, 0.9, 0.99. We note that E(T r ) and V(T r ) 
increase with p R for any fixed £. 

We conclude by pointing out that the purpose of the present note was to 
establish the quantitative foundations to a viable way to include refractoriness 
in neuronal diffusion models. Implementation of our approach to data analysis 
will be the object of future endeavors. 
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£ 


ti(S\ P ) 




E(T r ) 
Pr = 0-1 




E(T r ) 
P R = 0.5 




E(T r ) 
Pr = 0-9 




E(T r ) 
p R = 0.99 




0.5 


3.768002 Eh 


-2 


4.103229 


E+15 


3.692906 


E+16 


3.323615 


E+17 


3.655977 Eh 


-18 


1.0 


8.029989 Eh 


-1 


2.425535 


E+7 


2.182981 


E+8 


1.964683 


E+9 


2.161152 EH 


-10 


1.5 


4.661249 Eh 


-1 


4.020549 


E+4 


3.618494 


E+5 


3.256645 


E+6 


3.582309 EH 


-7 


2.0 


3.470051 Eh 


-1 


1.573636 


E+3 


1.416272 


E+4 


1.274645 


E+5 


1.402110 EH 


-6 


2.5 


2.866867 EH 


-1 


2.204449 


E+2 


1.984004 


E+3 


1.785603 


E+4 


1.964164 EH 


-5 


3.0 


2.502681 Eh 


-1 


5.871921 


E+l 


5.284729 


E+2 


4.756256 


E+3 


5.231882 Eh 


-4 


3.5 


2.258692 Eh 


-1 


2.264293 


E+l 


2.037863 


E+2 


1.834077 


E+3 


2.017485 Eh 


-4 


1.0 


2.083633 Eh 


-1 


1.102071 


E+l 


9.918635 


E+l 


8.926772 


E+2 


9.819449 EH 


-3 


4.5 


1.951789 Eh 


-1 


6.271095 


E+0 


5.643986 


E+l 


5.079587 


E+2 


5.587546 EH 


-3 


5.0 


1.848842 Eh 


-1 


3.983514 


E+0 


3.585162 


E+l 


3.226646 


E+2 


3.549311 EH 


-3 



Table 5: Feller model with ■& = 5, g = -70, v = -80 and £ = 0.5 • i (i = 1, 2, . . . , 10). 
In second column we have listed the FPT mean ti(S\p) with S = —50, whereas in 
columns three, four, five and six we have respectively listed the mean of refractoriness 
period for p R = 0.1, 0.5, 0.9, 0.99. 





V(S\ P ) 






V(T r ) 
Pr = 0.1 




V(T r ) 
Pr = 0-5 




V(T r ) 
Pr = 0-9 




V(T r ) 
p R = 0.99 


0.5 


1.395404 


EH 


-5 


1.683649 E+31 


1.363755 


E+33 


1.104642 


E+35 


1.336617 E+37 


1.0 


6.372482 


EH 


-3 


5.883257 Eh 


-14 


4.765411 


E+16 


3.859980 


E+18 


4.670576 E+20 


1.5 


2.241795 


EH 


-3 


1.620153 EH 


-9 


1.309681 


E+ll 


1.060603 


E+13 


1.283297 E+15 


2.0 


1.304116 


EH 


-3 


2.585121 Eh 


-6 


2.015619 


E+8 


1.625601 


E+10 


1.966008 E+12 


2.5 


9.313963 


EH 


-2 


6.143569 Eh 


-4 


4.051826 


E+6 


3.198777 


E+8 


3.85908 E+10 


3.0 


7.390905 


EH 


-2 


6.492913 Eh 


-3 


3.066862 


E+5 


2.286843 


E+7 


2.739949 E+9 


3.5 


6.238662 


Eh 


-2 


1.592475 Eh 


-3 


5.124507 


E+4 


3.451142 


E+6 


4.079672 E+8 


4.0 


5.478171 


Eh 


-2 


6.147816 Eh 


-2 


1.427634 


E+4 


8.366949 


E+5 


9.684436 E+7 


4.5 


4.940875 


eh 


-2 


3.064866 Eh 


-2 


5.588538 


E+3 


2.795397 


E+5 


3.144373 E+7 


5.0 


4.541290 


eh 


-2 


1.789263 Eh 


-2 


2.751622 


E+3 


1.172086 


E+5 


1.272925 E+7 



Table 6: For the Feller model and for the same choices of parameters of Tabled] in the 
second column we have listed the FPT variance V(S\q) with S = —50 and g = —70, 
whereas in columns three, four, five and six we have respectively listed the variance of 
refractoriness period for p R — 0.1, 0.5, 0.9, 0.99. 
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